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We present the result of three BeppoSAX observations of the X-ray binary pulsar A0535+26 in its quiescent 
state. The source is quite well detected up to 200 keV (6a detection in the stronger observation). No Iron line 
is detected in the MECS data (3a upper limit on its equivalent width of 150 eV). There is evidence of a soft 
excess below 2 keV. Pulsation is detected in all energy bands up to 10 keV, with a pulsed fraction of ~ 0.5, not 
varying with energy. There is a marginal detection (4a) of a cyclotron resonance feature (CRF) at 118 ±20 keV in 
the PDS spectrum. During the BeppoSAX observations HDE 245770, the optical counterpart to A0535+26, was 
monitored spectroscopically and photometrically. These observations show that the H a and Hp lines previously 
observed in absorption returned in emission, indicating that the Be disk formed again. The presence of pulsation 
and a CRF is a clear indication that accretion onto the polar caps is occurring, and that the propeller mechanism 
is leaking. The presence of a soft excess could be explained either in terms of thermal emission (kT ~ 1.2 keV) 
from the neutron star surface, as observed in other accreting X-ray pulsars, or in terms of an overabundance of 
Mg in the circumstellar matter. 



1. INTRODUCTION 

The class of High Mass X-ray binaries con- 
taining a Be star (Be/HMXRBs) is characterized 
by transient X-ray outbursts, sometimes reach- 
ing Eddington luminosities, interleaved by inter- 
vals of low luminosity (L x < 10 35 erg/s). It is 
possible to distinguish two types of outbursts pQ: 
Type I outbursts, lasting few days and due to the 
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increase of the accretion rate onto the neutron 
star as it approaches the periastron along the ec- 
centric (e > 0.5) orbit; and Type II outbursts 
(also called "giant" outbursts), lasting tens of 
days and thought to be due to sporadic episodes 
of mass ejection from the Be star. The details 
of the accretion onto the neutron star are com- 
plicated by the poorly known outflow mechanism 
in Be stars: because of the strong spinning of 
these stars, close to 70% of their break-up ve- 
locity 2 , circumstellar material will accumulate 
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in the equatorial region in the form of a quasi- 
Keplerian disk. Type I outbursts are thought to 
originate when the neutron star passes through 
this slow, dense disk. Quite recently, theoretical 
works |3I4| have shown that tidal interactions can 
produce truncation of the Be disk at a distance 
smaller than the periastron separation, suppress- 
ing the X-ray activity, and therefore explaining 
the missing of some Type I outbursts. In this 
model, long term (~ years) cycles arise as stellar 
material accumulate in the Be disk until tidal in- 
teractions sweep it away. The precession/ warping 
of the Be disk accounts for the irregularity of the 
Type II outburst recurrences. 

The study of the quiescent emission in 
Be/HMXRBs can give important clues on the ac- 
cretion flow geometry because of the tremendous 
difference in luminosity in the quiescent and out- 
burst states (more than three orders of magni- 
tude). In particular, we expect that the region 
responsible for the X-ray emission passes from 
being optically thick at high luminosity to be op- 
tically thin in quiescence, because of the forma- 
tion of a collision-less shock above the polar caps 
of the neutron star during outburst states 
These variations in X-ray luminosity are respon- 
sible of dramatic variations in the pulse shape , 
because for L x > 10 36 erg/s X-ray photons es- 
cape mainly in a direction perpendicular to the 
magnetic field lines (fan beam emission pattern), 
while at lower luminosities photons are emitted 
mainly in a direction parallel to the magnetic 
field lines (pencil beam). The comparison of pulse 
shapes at different luminosities (quiescent vs out- 
burst phase, but also at different orbital phases 
during quiescence) can give therefore some clues 
on the emission pattern geometry. 

Furthermore, the study of the quiescent state 
in Be/HMXRBs can give important information 
on the properties of mass accretion at very low 
accretion rates, and on the possible onset of a 
propeller state (see SectionQ]below), in which the 
accretion onto the neutron star is centrifugally 
inhibited. 

In this paper we present the results of a cam- 
paign of observations performed by BeppoSAX 
in order to study the quiescent emission of 
A0535+26 at different orbital phases, and there- 



fore at different luminosity levels. Simultaneous 
optical observations of HDE 245770 performed at 
the 1.5m Cassini Telescope of the Bologna Astro- 
nomical Observatory are also presented. 

2. A0535+26 

A0535+26 is the prototype of the 
Be/HMXRBs. It was discovered during a Type II 
outburst in 1975 |7I8| . and since then three other 
giant outbursts and about a dozen of normal 
outbursts have been observed UJ. X-ray pulsa- 
tions at ~104 s have been detected both during 
outburst ^01 an d in quiescence |1 lj . Because of 
the long database of observations it was possible 
to resolve the binary system ^2]> that exhibits 
an orbital period of 111 days, and an eccentricity 
of 0.47. The optical counterpart to A0535+26 is 
HDE 245770 [T3], of spectral type O9.7IIIe-B0Ve 
|14j . The Be star presents UBV photoelectric 
variability correlated to the X-ray behavior of 
A0535+26 [To\: in particular the V luminosity 
reached its maximum in 1975 (V ~ 8.8 m 
correspondence of a giant X-ray outburst. The 
minimum was observed in 1998 (V ~ 9.1 |lf)j). 
when the optical spectrum showed a dramatic 
change, with the emission lines slowly changing 
from emission to absorption. This phenomenon 
is interpreted in terms of Be disk loss J7| . 

3. OBSERVATIONS AND RESULTS 

3.1. X ray observations 

A0535+26 was observed three times at differ- 
ent orbital phases by the Narrow Field Instru- 
ments (NFIs) aboard the Italian-Dutch satellite 
BeppoSAX PH: LECS (0.5-10 keV [H), MECS 
(1.8-10 keV [23), and PDS (15-200 keV [21]). 
During the whole campaign the HPGSPC [221 
was switched off. The log of the observations, 
together with the observed NFIs count rates are 
reported in Tabled The first two observations, 
in which the neutron star approaches the peri- 
astron, show the same intensity, while the third 
observation, after the periastron passage, has an 
almost double count rate in all the NFI energy 
bands. 
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Table 1 

Log of the observations of A0535+26 performed by BeppoSAX 





OP 


Start Time 


Length 


(7) (10- 2 C/s) 




Orbital 




# 


(UT) 


(sec) 


LECS 
(0.5-2) 


MECS 
(2-10) 


PDS 
(15-200) 


Phase 


OBS1 


09736 
09738 


04/09/00 05:14:05 
05/09/00 00:41:56 


68865^ 
50020^ 


> 0.56 ±0.06 


3.2 ±0.1 


20 ±6 


0.606-0.618 


OBS2 


09923 
09924 

10835 


05/10/00 00:42:46 
06/10/00 04:29:26 

05/03/01 22:52:38 


100000" 
94584^ 

100000 ' 


> 0.56 ±0.05 


3.1 ±0.1 


11 ±5 


0.883-0.903 


OBS3 


10836 
10837 


07/03/01 02:39:28 
08/03/01 06:25:58 


100000 
40415 t 


> 1.20 ±0.05 


7.8 ±0.1 


24 ±4 


0.251-0.276 



a Ephemeris from Finger (private communication): T N = 45948.8(16) + N x 111.07(7) 




Figure 1. 500 s background subtracted MECS light curves of the third BeppoSAX observation in two 
energy bands (first two panels) and their ratio (lower panel). There is no appreciable spectral evolution 
during the whole observation 
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Beat Period: 103.40 +/- 0.02 s 




Start Time 11973 22:56:37:001 Stop Time 11076 17:38:16:606 

Figure 2. Epoch folding search result on the third 
A0535+26 observation. A Gaussian fit to the 
curve was used to compute the error on the mea- 
sured spin period 



3.1.1. Timing analysis 

During each BeppoSAX observation the source 
was quite stable, with time variability on time 
scales of the order of hundreds of seconds and in- 
tensity variability within a factor two (see Fig. ^ 
with the OBS3 light curve as observed in two 
MECS energy bands. Gaps in the data are due 
to passages in the South Atlantic Geomagnetic 
Anomaly, when the NFIs are switched off). The 
(4 — 10)/ (2 — 4) keV hardness ratio does not show 
any variation during the observations. 

After correcting photon arrival times to the so- 
lar system barycenter, MECS data were searched 
for periodic signals by performing an epoch fold- 
ing search and the result is shown in Fig. [21 for 
OBS3. Data were not corrected for the A0535+26 
orbital motion because of the short duration of 
the observation with respect to the orbital period 
of 111 days. Our result of 103.40 ± 0.02 s is in 
agreement with the other quiescent pulse period 
measurements |9I23) . confirming that the source 
is accreting material and channeling it onto the 
neutron star polar caps even when in quiescence. 

In Fig. we show the LECS and MECS pulse 
profiles as a function of energy for all the three 
observations (we were not able to extract pulse 
profiles from the PDS instrument, because of the 
very low statistics). In the first observation the 



2 4 6 8 10 

Energy (keV) 

Figure 4. Pulse fraction as a function of energy 
for OBS3. The dotted line represents the pulse 
fraction measured on the whole 4-10 keV energy 
range 



source displays a double peak profile, typical of 
the high energy band 24 6 . In OBS2 the pulse is 
clearly single peaked, of the same kind observed 
by RXTE at orbital phase 0.64 [73 , while in the 
third it is structured, with the presence of a pro- 
nounced first peak and a smaller second one. The 
pulse fraction within the uncertainties, does not 
depend on energy, as it is shown in Fig. 01 for 
OBS3. Our result is in agreement with the pulse 
fraction of 0.53 measured by |23| . 

3.1.2. Spectral analysis 

The LECS and MECS energy spectra for all 
the observations were extracted using standard 
procedures, as outlined in |25j . The combined 
LECS/MECS spectra are well fit both by an ab- 
sorbed power law and a thermal bremsstrahlung. 
Also a blackbody fit is statistically acceptable, 
but the resulting emitting region is unphysically 
small (see on this topic), and no absorption 
is required. These results are in agreement with 
and are listed in Tabled In Fig. we show 
the power law fit for the third observation, in 
which a soft excess below 2 keV is quite evident. 
This excess can be modeled both with the addi- 
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Figure 3. Background subtracted LECS and MECS pulse profiles for all the three A0535+26 observations 
in four energy bands: 0.5-2 keV (LECS; first row), 2-4 and 4-10 keV (MECS; second and third rows). 
In the last row it is shown the 2-10 keV pulse profile. All the profiles are phased at MJD 11791 
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Table 2 

A0535+26 LECS+MECS spectral fits' 



0BS1 OBS2 OBS3 



Power Law 

N ff (10 22 cm 
a 

Xl (dof) 



1.10 ±0.19 
1.98 ±0.11 
0.90(68) 



1.03 ±0.13 
1.92 ±0.09 
1.17(98) 



1.33 ±0.09 
1.85 ±0.05 
1.33(208) 



Bremsstrahlung 

N H (10 22 cm- 2 ) 
kT (keV) 
Xl (dof) 


0.75 ±0.14 
7.4 ± 1.2 
0.90(68) 


0.75 ±0.10 
7.7 ± 1.1 
1.07(98) 


1.02 ±0.07 
9.4 ±0.8 
1.17(208) 


Black Body 

kT (keV) 
R (Km) 
Xl (dof) 


1.09 ±0.03 
0.088 ± 0.004 
1.37(69) 


1.14 ±0.02 
0.080 ±0.003 
1.23(99) 


1.26 ±0.01 
0.106 ±0.002 
1.30(209) 


L x (0.5-2 keV) b 
L x (2-10 keV) 6 


0.6 
1.5 


0.5 
1.5 


1.2 
4.4 



a Uncertainties at 90% confidence level for a single parameter 

b Un-absorbed luminosities in 10 33 erg/s assuming a distance d = 2 kpc [2*7| 



model: wabs pow 



BeppoSAX OBS3 Observation of A0535+26 




Energy (keV) 

Figure 5. BeppoSAX A0535±26 count rate spectrum (plus signs; LECS in black and MECS in gray) and 
power law best-fit continuum (histogram), together with the fit residuals. Fitting parameters are listed 
in Tabled 
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Table 3 



model: pow gausaba 



BeppoSAX 0BS1 Observation of A0535+26 



Figure 6. PDS 15-200 keV count rate spectrum 
for OBS1. Top panel: count rate spectrum (plus 
signs) and best fit continuum model (histogram) . 
Bottom panel: residuals to the power law plus 
Gaussian in absorption model. The CRF nor- 
malization was put to zero, in order to show the 
feature significance 



tion of a black body component (kT — 1.2 keV; 
X 2 /dof = 213.5/206), or by an over-abundance of 
Mg in the circumstellar matter (14 ± 3 times the 
Solar abundance; x 2 /dof = 232.5/208). Our data 
do not show the presence of an iron iC-shell line in 
6.4-6.9 keV: the 3<r upper limit on its equivalent 
width is 150 eV. 

For extraction of the PDS net spectra we used 
the standard procedure [2S1 f° r a U the observa- 
tions but the third where one of the offset fields 
used for the determination of the background was 
contaminated by the Crab nebula, and therefore 
could not be used in the pipeline. In all the ob- 
servations A0535+26 is clearly detected in the 
PDS (see Table QJ up to 200 keV. To our knowl- 
edge, this is the first detection up to these energies 
of the quiescent state of A0535+26. Because of 
the low statistics only a power law model could 
reasonably be applied for fitting the PDS spec- 
tra. The best fitting results are listed in Table |3j 
These results indicate a hardening of the spec- 
trum for E>15 keV consistent with a = 1.1 — 1.4 
observed by HEXE (15-200 keV) on-board MIR 
during an outburst UJ. 



A0535+26 PDS 


power law 


spectral fits 






OBS1 


OBS2 


OBS3 


a 


1.0 ±0.5 


2.1 ±0.8 


1.0 ±0.2 


xl (dof) 


0.55(8) 


0.86(8) 


0.86(8) 



CRF best fit parameters 



£ C rf (keV) 
a (keV) 
E.W. (keV) 



5 +0 - 4 

U -°-0.8 
+21 
-19 
20 



118 

20_ 19 
811S 



The PDS spectra for OBS2 and OBS3 do not 
show any deviation from the power law fit, while 
the OBS1 spectrum presents structured residu- 
als in the range 100-200 keV (see the slightly 
worse xt f° r OBS1 in Table 0). Because of a 
previous observation by OSSE of a cyclotron 
resonance feature (CRF) at 110 keV (observed 
during a giant outburst), we tried to model the 
feature with a Gaussian in absorption |30|. as 
it has been done with other CRFs observed by 
BeppoSAX. We prefer this kind of model with 
respect to the Lorenzian model |31| because the 
latter has the disadvantage of modifying the con- 
tinuum (see discussion in [6'2\). In Fig. we show 
the power law plus Gaussian in absorption fit for 
OBS1, where in the residuals panel we put the 
CRF normalization to zero in order to show its 
significance. The CRF best fit parameters are 
listed in Table 03 in agreement with the OSSE 
results. 

In order to check whether the observed spectral 
feature was real or due to an instrumental effect, 
we performed a normalized Crab ratio |33j on the 
PDS count rate spectrum of the first observation. 
This technique has been successfully applied for 
the detection of CRFs in other X-ray pulsars [3"2~] 
and consists of the ratio between the source count 
rate spectrum and the count rate spectrum of the 
Crab Nebula. As this second spectrum is known, 
with great accuracy, to be free of features and to 
be modeled at first order with a power law in a 
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A0535+26 0BS1 Crab Ratio 
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Figure 7. Normalized Crab ratio for the first PDS 
observation of A0535+26. First panel: ratio be- 
tween the A0535+26 and Crab count rate spectra. 
Middle panel: Crab ratio multiplied by the func- 
tional form of the Crab spectrum, a E -2,1 power 
law. Bottom panel: The result from the previous 
operation divided by the functional describing the 
continuum (as listed in Table |3J) 



very broad energy range, this ratio is quite well 
suited to enhance the presence of features in the 
spectrum. Furthermore the ratio is in first ap- 
proximation independent from the calibration of 
the instrument. In order to enhance the devi- 
ations from the continuum we multiply the ra- 
tio by a E -21 power law, that is the functional 
form of the Crab spectrum, and we divide by the 
functional describing the continuum shape of the 
source (from this the name normalized Crab ra- 
tio). The procedure is described in Fig. where 
we plot the result of each different step used to 
obtain the final result in the case of A0535+26. 

3.2. Optical observations 

The optical counterpart to A0535+26, the Be 
star HDE 245770, was monitored at the 1.5m 
Cassini telescope of the Bologna Astronomical 
Observatory almost continuously since its discov- 
ery in 1975 |33 . During the BeppoSAX cam- 
paign several medium-resolution optical spectra, 



with dispersion ranging from 1.8 to 4 A/pixel and 
covering the 4000-8500 A range, were acquired 
simultaneously with the X-ray pointings by us- 
ing the BFOSC instrument. One of these spectra 
is shown in Fig. [SJ In all of them it is clearly 
seen that H a and H^ are in emission, that is the 
Be star spectrum returned to its normal Be spec- 
trum. This is in agreement with the rebuilding of 
the Be circumstellar disk around HDE 245770. 

At the time of these observations, optical pho- 
tometry data collected at the nearby Loiano 60 
cm telescope shows that the star was at magni- 
tude v ~ 9.0 [351- 

4. DISCUSSION 

Recent RXTE observations [23] of the quies- 
cent emission from A0535+26 at L x ~ 4 x 10 33 
erg/s (2-10 keV; two orders of magnitude lower 
than that observed in past quiescent states |11|1 
pose serious problems for the standard accretion 
theory (see, e.g. [3S])- Indeed, at this luminos- 
ity level the corotation radius (the distance from 
the neutron star at which the centrifugal force 
equals the gravitational pull) is smaller than the 
magnetospheric radius (the distance at which the 
magnetic pressure equals the ram pressure), and 
therefore matter cannot enter the magnetosphere 
and being accreted because of centrifugal forces 
(this is the so called propeller state (33)- But the 
detection of X-ray pulsations clearly contradicts 
that the source is in a centrifugal inhibited state. 
Optical observations before and during the RXTE 
observations also showed that the JHK infrared 
excess and IR/optical lines emission were missing 
|17j . a clear indication of a Be disk loss episode 
in HDE 245770. 

In order to accommodate this discrepancy 
Ikhsanov j3EJ proposed that the mechanism of 
plasma penetration into the magnetosphere is 
not through standard interchange instabilities, 
but through magnetic reconnection of the field 
lines. In this way, both the presence of pulsation 
and the small polar cap radius derived from the 
RXTE spectral fits are explained. 

Our X-ray observations confirm the low lumi- 
nosity state of the source, and therefore its cen- 
trifugally inhibited state. The LECS observations 
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Figure 8. Optical spectrum of HDE 245770 taken at the 1.5m Cassini telescope of the Bologna Astro- 
nomical Observatory. The star returned to its normal Be spectrum 



were able to clearly detect a soft excess that could 
be explained both in terms of an additional ther- 
mal component, as already observed in other ac- 
creting X-ray binary pulsars, like Her X-l |39) . 
or overabundance of Mg (~ 10 times the Solar 
abundance) in the proximity of the source. The 
almost absence of spectral variability during the 
orbit seems to point to the absence of any (tem- 
porary) accretion disk, as observed during a gi- 
ant outburst phase from the detection of spin and 
QPO frequency variations |12j . 

The main difference between BeppoSAX and 
RXTE observations is that the Be disk around 
HDE 245770 returned in place. We would have 
expected a higher X-ray luminosity with respect 
to the disk-less observations, due to the greater 
amount of circumstellar matter that can be ac- 
creted. Therefore the presence of the Be disk 
seems not to alter the accretion rate onto the neu- 
tron star. 

The high energy observations indicate a spec- 
tral hardening, which is in disagreement with the 



X-ray spectra observed during the A0535+26 ac- 
tive phase, characterized by a spectral softening 
above 25 keV |29|. This is further evidence that 
the kind of accretion in the quiescent and the ac- 
tive phase seems different. In the first BeppoSAX 
observation there is a marginal (4a) evidence of 
a CRF at 110 keV, but not of a 55 keV feature 
(reported by |40|). This is a confirmation of the 
OSSE observation of the CRF detected during 
the 1994 giant outburst. The same CRF energy 
observed in two completely different luminosity 
states constrains the physical conditions in the 
region in which the CRF originates to not de- 
pend on the accretion rate, and this poses serious 
limits on the actual emission models. 
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